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> : 

O ■ Abstract: We present an update of the JetViP 1.1 program for performing 
^ ■ fixed NLO calculations in jet production including direct and resolved compo- 
rt '. nents in a continuous range of photon virtuality Q^. The new version allows 

I to access the full event record on the parton level. The program is set up 

>- I such that hbook can be used to fill histograms. The phase space generator has 

^ ■ been optimized and the azimuthal dependence of the cross sections is taken 

. into account in LO. We comment on recent comparisons between various NLO 

^ ! programs for jet production at HERA. We demonstrate that the J2 Et cut for 

I dijet cross sections is not infrared safe. 
0^ ■ 

^: 

1 Introduction 

D , 

"T^ ! Several fixed order programs are on the market now for calculating jet cross sections for 

; non-zero photon virtuality at HERA, as there are D I SENT DISASTER++ [|], MEPJET 

^ ■ and JetViP which is based on the calculations in ||, |[. The first two programs are 

■ based on the subtraction method, whereas the latter two employ the phase-space-slicing 
method. In the photoproduction limit 0, the deep-inelastic scattering (DIS) matrix 

elements with a direct coupling of the photon to the partons from the proton show an 
additional initial state singularity on the photon leg, which has to be subtracted and 
absorbed in the photon structure function. Furthermore, one has to calculate a so-called 
resolved component, where the photon serves as a source of partons. These direct and 
resolved photon-proton scattering processes have been calculated by three groups in NLO 
QCD 0, H, ||. The calculation of uses the subtraction method, whereas the other two 
groups 1^, use the phase-space slicing method. One of the main features of JetViP is 
the possibility to include a resolved virtual photon component in NLO. In this way the 
photoproduction limit can be taken, which allows to cover the full range in accessible 
at HERA. 

The phase-space generator of JetViP was originally built around the i^^ep-niodified cone- 
algorithm. Though the phase-space could be optimized for this specific case, the program 
was obviously limited in its applicability and furthermore it was rather inconvenient to 
calculate spectra of variables other than the jet transverse energy and rapidity. Therefore, 
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in the new version partons rather than jets are accessible and an event record is provided 
in a user-routine, i.e., the four- vectors of all particles involved in the process are available 
together with the basic kinematical variables such as Q^, xbj, y, ■ ■ ■ and the weight from 
the matrix elements convoluted with the parton distribution functions. The analysis of 
the events is then left entirely to the user routine. 

Two further important improvements have been achieved. First, the phase space generator 
has been rewritten, as to generate events in {y, Q^)-space, which greatly improves the 
statistics. Second, the azimuthal dependence of the matrix elements, which has some effect 
for cuts in the laboratory- frame, have been implemented in LO. As we will demonstrate, 
this is sufficient for small and when the NLO corrections are not too large. 

Recently, a comparison of the above mentioned NLO programs for DIS has become avail- 
able flOl- We have recalculated the cross sections produced by JetViP presented in [|l^ 



with the new version. As will become clear, differences in the cross sections produced 
by the programs are due either to statistical problems or to infrared (IR) sensitive cuts 
on dijet variables. Similarly, a comparison of photoproduction calculations has been per- 
formed [|ll[|. We compare photoproduction cross sections from JetViP to those from one 



of the programs used in |TT| and briefly comment on the implications for the DIS region. 



Before discussing these more physical topics, we describe the new user-routine with 
common-blocks and the new steering file. Further details of the JetViP program can 
be found 11. 



2 Program changes 
2.1 user routine 

The events generated by JetViP have to be analyzed in a routine called user, 
subroutine user (if lag, weight , if ill) 

real*8 vgweight , weight , Qs , shad , Ws , xbj , y , xa , xb , plab , pcms 
integer if lag, if ill 
common/parton/plab (7,7), pcms (7,7) 
common/kin/Qs , shad , Ws , xbj , y , xa , xb 
common/vgplot/vgweight 



• if lag: the routine user is called once at the very beginning of JetViP for initializa- 
tion with if lag=l and weight=OdO. After this initialization, events are generated 
with if lag=2. At the very end user is called for termination with if lag=3 and 
weight=OdO. 
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• weight: weight of the event, i.e., the matrix elements convoluted with the PDF's 

from the proton and/or photon, together with running Uqed and (in pb). The 
entry for the histogram has to be further multiplied with the weight from VEGAS, 
which is vgweight. The histogram entry therefore has to be vgweight*weight. 

• if ill: for if ill=0 the phase space is scanned by VEGAS with low statistics. In this 
stage the histograms must not be filled. For if ill=l the histograms can be filled. 

• plab(7,7): matrix containing the four- vectors, Et, rj and for 7 different particles 

in the laboratory frame. The 2;-dircction is given by the incoming proton or, for the 
7*7 mode, by the real photon. 1st entry: particle. 2nd entry: variable. 



/ plabd.i) \ 




^ incoming electron ^ 




/ plab(j,l) \ 




( \ 


plab(2,i) 




outgoing electron 




plab(j,2) 




7] 


plab(3, i) 




incoming proton 




plabCj ,3) 







plab(4,i) 




incoming photon 




plab(j,4) 




Px 


plab(5,i) 




outgoing parton 1 




plabCj ,5) 




Py 


plab(6,i) 




outgoing parton 2 




plabCj ,6) 




Pz 


plab(7,i) y 




outgoing parton 3 j 




plabCj ,7) ) 




[ E ) 



The outgoing partons are not ordered in Et- The incoming photon vector, i.e., 
plab(4,i), is set equal to zero. 



• pcms(7,7): matrix containing the four- vectors, EtiTi and for 7 different particles 
in the hadronic, i.e., 7*P cms frame. In the 7*7 mode the 7*6 cms frame is chosen. 
The ^-direction is given by the incoming proton/real photon. Same conventions as 
for plab(7,7). The incoming and outgoing electron vectors, i.e., pcms(l,i) and 
pcms(2,i), are set equal to zero. 

• Qs,shad,Ws: photon virtuality = — g^, total cms energy sh = 4:EaEb and 
hadronic cms energy W'^ — {P + q)"^. 

• xbj,y,xa,xb G {0,1}: Bjorkcn-a; variable Xbj, y-variable, momentum fraction of 
the parton in resolved photon Xa and momentum fraction of the parton in proton 

Xb- 

2.2 Steering file 

An example of the steering file is given in appendix B. 
Contribution 

• iproc: integer e {1, 2}. Type of process. 1 = eP, 2 = e+e~ (7*7-mode) 
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• isdr: integer e {1,2,3,4}. Selection of the D, SR, SR* or DR contribution. 1 =D; 
2 =SR; 3 =SR*; 4 =DR. Only one of the contributions D, SR, SR* or DR can be 
calculated at a time. 

• iborn: integer e {0, 1}. Born process (2 —>■ 2). 

• itwo: integer e {0, 1}. NLO twobody corrections (2^2). 

• ithree: integer e {0, 1}. The 2 — > 3 contributions above the i/g cut. 

• isplit: integer e {0, 1}. Initial state singularity for the virtual photon (2^2). 

• iqcut: integer G {0, 1}. Insert a finite cut-off t/s (cutmin) into the integration on 
the virtual photon side for the 2 — > 3 NLO contributions (= 1) or not (=0). Setting 
iqcut=l is necessary for photoproduction. 

In the following the main settings for LO and NLO cross sections are summarized in 
tabular form. These settings hold for all 4 values of isdr. 



Initial State 

• Ea, Eb: real*8. Energies of the incoming lepton a and hadron/lepton b in GeV. 

• iwwa: integer e {0, 1}. Selects formula for the photon flux in the case of photopro- 
duction. For iwwa=0 the cross section is intergated analytically up to Q'^ax^ whereas 
for iwwa=l it is integrated using the value of thmax below. 

• thmax: real*8. Maximum scattering angle of electron in photoproduction for iwwa=l. 

• P2max: real*8 > 0. Maximum value of (7*7 mode). 

• iwwb: integer e {0, 1}. see iwwa (7*7 mode). 

• thetbmax: real*8. Maximum scattering angle of electron for iwwb=l (7*7 mode). 



LO NLO (PIS) NLO (Q^ = 0) 7* ^ qq 



iborn 

itwo 

ithree 

isplit 

iqcut 



11 10 

1 10 

1 10 

1 1 

10 
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Subprocess 



Nf: real*8 G {1,2,3,4,5}. Number of active flavours, 
lambda: real*8 > in GeV. Value of Aqcd for Nf flavours. 

ialphas: integer G {1,2}. One-loop (ialphas= 1) or two-loop (ialphas= 2) 
formula for the strong coupling constant as without thresholds. For ialphas = 3 
the value of is taken from the PDFLIB (automatically adjusts A). 

ycut: real*8 > 0. Value of i/s- The independence of the NLO cross sections on ys 
has been tested for yg G [10~^, 10~^]. Large statistical errors occur for too small ys- 
Recommended value: ys = 0.5 • 10"'^. 

idisga: integerG {0, 1}. Since all partonic cross sections in JetViP are imple- 
mented in the MS-scheme, a photon PDF constructed in the DIS^ scheme has to be 
transformed into the MS-scheme. This is done for the photon on side a by setting 
idisga= 1 [|]. 

ipdftyp: integer G {1,2,3,4}. Selects the PDF for the resolved virtual photon. 
1 = SaS PI, 2 = GRS 111, 3 = PDFLIB ^ (for real photons with = 0). 



igroupa: integer. For ipdftyp= 3,4 this selects the group from the PDFLIB (see 



manual [T^). For ipdftyp= 1 this represents isassetG {1,2,3,4}, which selects 



input scale and scheme of SaS virtual photon PDF (see |T! 
iseta: integer. For ipdftyp= 3,4 this selects the set from the PDFLIB (see manual 



1^). For ipdftyp= 1 this represents isasp2G {1,2,3,4,5,6,7}, which selects the 
scheme used to evaluate off-shell anomalous component in SaS virtual photon PDF 
(see 



II)- 



• igroupb: integer. Selects the group from PDFLIB for the resolved component for 
particle b (see manual [ll4|). 

• isetb: integer. Selects the set from PDFLIB (see manual |T^ ). 

• a, b, c: real*8 > 0. Choosing the overall-scale fi'^ according to /i^ = a + bQ'^ + cp'^. 
The pt is the largest E^^ of an event. 

Phase space integration 

The phase space of JetViP is generated in the following variables. 

• Q2min, Q2max: real*8 > 0. Minimum and maximum value of the photo virtuality 
Q^. Setting Q2niin =0.dO will produce a photoproduction cross section. 

• ymin, ymax: real*8 G {0, 1}. Minimum and maximum value of y. 

• ptmin: real*8 > 0. Minimum value of Et on parton level. Set to ^E^?, 



1 ipjet 

min' 
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Vegas and Output 

• ipoint: integer > 0. Defines the number of events generated in tlie pliase space 
defined above. Typical value: ipoint = 10.000.000. In NLO for smaller Us larger 
compensations occur and therefore the statistical errors become larger. 

• itt: integer > 0. Number of iterations for Vegas. Recommended value: itt= 5. 

• jfileout: character*20. Name of the output-file in which the histograms are stored. 



3 Comparison with other programs 



3.1 Infrared safe dijet cross sections 



Recently, a comparison of NLO programs for the DIS region has become available | TI| , 
which shows a good agreement between the programs DISENT and DISASTER, systematic 
deviations of MEPJET, being typically of 5-8% lower than the other programs and finally 
the overall results of JetViP being comparable to DISENT and DISASTER. However, in 
some cases JetViP deviates from these programs by up to 20%. Furthermore, a significant 
dependence of the cross sections produced by JetViP on the slicing parameter Us has been 
reported. In the following, these results are discussed in detail and the numbers for most 



scenarios described in iTOl are recalculated. 



The details of the technical settings for the comparison can be found in . Basically, 
dijet cross sections have been calculated under HERA conditions in the Breit frame, using 



the inclusive /c^ algorithm ||T^. The authors define as a central scenario the cuts 

30 < < 40 GeV^ 0.2 < y < 0.6, E]^'^ > 5 GeV (1) 

where and E^ are the highest and second highest Et of an event. Then a number of 
additional cuts are proposed in six scenarios, which are meant to ensure IR safeness of the 
jet cross sections. In the following we will concentrate on the comparison for scenarios 

1- 4, since scenarios 5 and 6 are not problematic and sufficient agreement between the 
different NLO programs has been found. As we will show later, the scenario 1 contains IR 
sensitive cuts. In this section we concentrate on the IR safe scenarios, namely scenarios 

2- 4, which are collected for reasons of completeness in the Tab. [I|. 

Recalculating the scenarios 2-4 with the new JetViP version 2.1 using higher statistics and 
the optimized phase space generator we find the results shown in appendix A, which mostly 
confirm the numbers from [|10]. Some numbers however show better agreement. JetViP 
agrees with DISENT and DISASTER on the 1% level within the statistical errors. Only 
for scenario 2d JetViP lies 2% below the D I SENT/DISASTER results. This is a statistical 
problem of the rather extreme scenario 2d. The cut on the higher Et jet is at 40 GeV, 
however all the E^^s down to 2.5 GeV have to be generated on the parton level, since they 
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No. 


^Tlmin/GeV 


No. 




QLx/GeV^ 


No. 




Z/max 


2 a) 


8 


3 a) 


3 


4 


4 a) 


0.01 


0.05 


2 b) 


15 


3 b) 


30 


40 


4 b) 


0.2 


0.6 


2 c) 


25 


3 c) 


300 


400 


4 c) 


0.9 


0.95 


2 d) 


40 


3 d) 


3000 


4000 









Table 1: Different IR safe dijet scenarios. The y and regions refer to an asymmetric 
cutting scenario with i^Timin = S GeV 



can contribute to the cross section after the jet recombination. Due to the steep fall-off 
of the Et spectrum, too few events are generated at large Erp- Overall, the agreement 
found between the three programs DISENT, DISASTER and JetViP is very good, which is 
especially promising since the programs are based on different methods to treat soft and 
collinear singularities. 

A reason for the systematic deviation of MEPJET from the other three programs, being 
overall about 5-8% too small, might be due to not partially fractioned matrix elements 
p!6| . In principle the partial fractioning is not necessary if one goes to extremely small 
Smin values. This, however, leads to large statistical errors and the calculation of precise 
results becomes impractical. Problems of this kind in connection with the phase space 
slicing technique are well known already from older works on jet production from e+e" 
anihilation (see e.g. 0]). 

A last point that has to be sorted out is the observed ys dependence of the JetViP cross 
sections. Except for the scenarios lb and Ic, which are discussed below to be IR sensitive, 
these dependences are a statistical problem. We have observed that the adaptation proce- 
dure implemented in the Monte- Carlo integration routine VEGAS becomes unreliable 
in the region of small slicing parameters ys and the errors are drastically underestimated. 
These kind of problems have been noticed recently and an alternative integration algo- 
rithm has been proposed [0. To avoid such problems we have switched off the adap- 
tation in VEGAS and simply generated an extremely large number of events. With this 
setting we have recalculated the ys dependence for three values of the slicing parameter. 



Vs 


scenario 2a 


scenario 3d 


10-3 


121.5±0.8 


1.99 ±0.02 


10-^ 


121.4±2.2 


2.0 ±0.05 


10-5 


116.1±6.5 


1.80 ±0.13 


DISENT 


119.5±0.3 


1.985±0.003 


DISASTER 


119.8±0.4 


1.998±0.003 



Table 2: The dependence of JetViP cross sections on the phase space slicing parameter 
ys for two scenarios. 
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Us = 10 ^, 10 ^ and 10 ^ for scenarios 2a and 3d which are among the critical cases in 
T0| . The results are shown in Tab. ^. No dependence on i/s within the statistical errors 



is observed. Of course, the errors become larger, the smaller i/s gets and it is very hard to 
calculate the cross section for values of i/s < 10~^ since the compensation between large 
positive and negative contributions lead to large statistical erros, which is a well-known 
problem of the phase-space slicing method. However, one can obtain reliable results with 
JetViP for values of i/s between 10~^ and 10~^ and in this region also use the adaptation 
technique implemented in VEGAS. 



3.2 Infrared sensitive scenarios 



The authors in |]T0| have proposed several ways to avoid IR sensitive regions for the dijet 
cross sections. These are the scenarios la-c, collected in Tab. R|. 



No. 


additional jet cut 


1 a) 


Erimm > 8 GeV 


1 b) 


Mjj > 25 GeV 


Ic) 


{Eti + Et2)>17 GeV 



Table 3: Proposed ways to avoid IR sensitive regions. 

Unfortunately, the scenarios lb and Ic are not IR safe. This can be understood in terms 
of the E^, E^ plane, shown in Fig. 1. The thick solid lines show the symmetric cuts 

5 GeV 



Ej^ > 5 GeV, which are known to be IR sensitive, due to the point E^ 
1^, The (negative) twobody contribution, containing the analytic NLO corrections 



Ej^ 



can not be compensated by (positive) threebody contributions due to the limited phase 



5 GeV— 







A cut / 




\ ' X > 



5 GeV EE Tout 



Figure 1: The available phase space for the two highest E^ jets. 
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76.4 


79.6 




0.250 




76.6 


79.7 




0.375 




76.1 


79.2 




0.500 




75.6 


78.1 




0.625 




73.7 


77.1 




0.750 




72.3 


75.5 




2.000 




54.2 


55.2 





Table 4: T/ie J2 Et scenario with an additional cut Ej-^ = 8.5 GeV+ A 

space for the third particle. To avoid this region, one can introduce a A-cut on the highest 
Et jet, E^ = ET,min + ^ with A ^ 1 GeV. Furthermore, the proposed scenario Ic is shown 
in Fig. 1. However, the point -E^ = E^ = 8.5 GeV corresponding to a symmetric scenario 
is not avoided by the J2 Et cut and not enough threebody phase space is added to this 
point. 

This IR sensitivity is the origin of the i/s dependence of JetViP. Furthermore, the cross 
sections obtained with two different methods, namely the phase space slicing method 
and the subtraction method, do not agree. Apart from the i/s dependence of JetViP, an 
explicit way to see the IR sensitivity is to introduce a A-cut in addition to the J2 Et 
cut, i.e., to demand an additional cut on the higher Et jet with E^ = ET,min + A GeV, 
where ET,min = 8.5 GeV. We have introduced such a cut and calculated the cross sections 
for various A values for the two programs DISENT and JetViP. The results are shown in 
Tab. ^. Physically meaningful cross sections will drop for rising A. However, as can be 
seen in Tab. ^ the NLO cross sections for the J2Et scenario still rise up to A = 0.375 
and only reach the original value of the cross section around A = 0.625. This behaviour 
is seen for both programs. We have checked several values of A and found A = 2 GeV 
sufficient to avoid IR sensitivity. As can be seen from the table, at this value the programs 
JetViP and DISENT agree at the 2% level. 

The same is true for scenario lb, where Mjj > 25 GeV. This cut washes out the IR 
sensitive region, but does not completely avoid it. This can again be seen in the ys 
dependence of the cross sections from JetViP. By introducing an additional A-cut, as for 
the J2 Et case, one would again see a rise of the NLO cross section with increasing A. 



3.3 Photoproduction 

Recently, a comparison of calculations in jet photoproduction has been performed [|n 
in which the authors find that the three calculations 0, H, Q agree within the statisti- 
cal accuracy of the Monte Carlo integration. Only in certain restricted regions of phase 
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space, the calculations differ by up to 5%. Since the program JetViP allows to perform 
the limit starting from the DIS matrix elements, a comparison with one of the 

photoproduction calculations will serve as a further consistency check of the two methods 
for extracting coUinear and soft phase-space regions. We have therefore calculated inclu- 
sive single-jet cross sections for photoproduction with JetViP and compared them to the 
calculations in [^]. The results are shown in Fig. ^ for the Et and rj spectra of the direct 
(dashed) and resolved (full) components in NLO. The agreement with the predictions 
from 1^ (full dots) is perfect for both components with differences below 1%. All other 
comparisons performed so far between the two programs show a similar good agreement. 



10' 



^ 10 

> E 



llO^r 

% 10 

■a E 



Res 
Dil- 



10 



12 



14 



16 




Figure 2: Single-jet inclusive direct (lower curves) and resolved (upper curves) jet cross 
sections with JetViP compared to the calculations /j^ (dots) for photoproduction. 



Since the matrix elements for the resolved component are taken from 0, it is not surprising 
that it agrees for both programs. The direct component in JetViP is, however, completely 
independent from the calculations in |p. The excellent agreement of the two programs 
therefore establishes the consistency between the calculations for photoproduction and for 
the DIS case for the phase space slicing method. Since the equivalence of the subtraction 
and the phase space slicing method has been shown in , we conclude that JetViP also 



agrees with the calculations based on the subtraction method in the photoproduction 
regime 0. 



3.4 Dependence on azimuthal angle 

We have implemented the full azimuthal dependence of the jet cross sections in the new 
JetViP program in LO according to the formula 

^ = v4 + 5cos</) + Ccos2<^ , (2) 
d(p 
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where denotes the azimuthal angle of the jets around the virtual photon direction in 
the hadronic cms frame. The lepton plane defines = 0. The coefficients A, B, C are 
related to the polarization density elements of the virtual photon. 

The dependence in @ can be integrated out if no cuts on the jets are imposed in the 
laboratory frame. This is, however, seldom the case experimentally. In the following we 
check how large the effect of using the integrated and unintegrated matrix elements is in 
LO. We calculate cross sections in LO within the central scenario defined above with no 
additional cuts and probe the region G [0, 100] GeV^. In Fig. |] three rapidity regions 
are shown in the laboratory frame, namely 

- 1 < r]iab < 0.5, 0.5 < r]iab < 1-5, 1.5 < r]iab < 2.8, (3) 

corresponding to the backward, central and forward region of the detector. We have 
plotted the ratio of the integrated over the unintegrated matrix elements, denoted as 
R. As one sees, for < 10 GeV^ the 0-dependence is below 2.5% for all regions and 
therefore negligible. In the central region the dependence stays below 5%, whereas in 
the more extreme regions the dependence can reach 10% for the largest Q^-bin shown. 



This is in accordance with the results in [211. We have further checked that the difference 
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Figure 3: The ratio of unintegrated over integrated LO cross sections in four different 
rapidity regions in the lab frame. 
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between integrated and unintegrated has a maximum at = 100 GeV^. For larger 
virtualities, the difference again diminuishes and falls below 5% for all rapidity regions 
above = 500 GeV^. If the NLO corrections stay reasonably small, e.g., around 20%, 
one can estimate the error due to cuts in the laboratory-frame in NLO to be below 2%, 
for any value of Q^. 

We conclude that for the main region of applicability of JetViP, namely the region of 

< 10 GeV^, where also the resolved cross sections plays a role, the 0-dependence 
is described correctly in the present JetViP version. For all other Q^, the dominant 
contribution has been taken into account and a correct description is given if the NLO 
corrections are not too large. In any case, the corrections will not be larger than 10%. 

4 Summary 

We have presented an updated version of the JetViP 1.1 program for jet production in 
eP- and 7*7-scattering. The most important changes are: 

• a full event record is provided on the parton level 

• the complete analysis of the events has to be done in the user routine 

• the phase space is generated in {y, (5^)-space 

• the (^-dependence has been implemented in LO 

We have numerically shown the equivalence of two methods used to extract singular phase 
space regions, namely, the subtraction method and the phase-space slicing method, by 
comparing 4 existing NLO programs. The programs DISENT, DISASTER and JetViP agree 
on the 1% level within the statistical errors of integration. Any differences observed 
before in these programs are either of statistical nature or due to IR sensitive cutting 
scenarios. In particular, we have demonstrated the IR sensitivity of the J2 Et scenario. 
The numerical equivalence of phase space slicing and subtraction method have furthermore 
been underpinned by comparisons in the photoproduction regime. 

The difference between using azimuthal dependent matrix elements in LO and those 
where this behaviour has been integrated out is negligible for all rapidity regions for 

< 10 GeV^. By taking into account the azimuthal dependence of the matrix elements 
in LO also cross sections involving cuts in the laboratory frame are described correctly in 
the present JetViP version for all practical purposes. 
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A Numerical results 



Results of the comparison between the programs DISENT , DISASTER , MEP JET and JetViP 
in the DIS region for IR safe scenarios 2-4. The numbers are taken from [0, except for 
the NLO JetViP numbers, which have been recalculated. 



scenario 


DISASTER++ 


DISENT 


JETVIP 


MEP JET 


2 a) 


119.8 ±0.4 


119.5 ±0.3 


121.5 ±0.8 


113.5 ±0.2 


LO: 


41.66 ibO.08 


41.77 ±0.06 


41.75 ±0.03 


41.72 ±0.03 


2 b) 


16.58 ±0.09 


16.53 ±0.05 


16.21 ±0.07 


15.74 ±0.08 


LO: 


6.19 ±0.02 


6.22 ±0.01 


6.214 ±0.005 


6.221 ±0.003 


2 c) 


2.08 ±0.03 


2.052 ±0.008 


2.04 ±0.1 


1.908 ±0.008 


LO: 


1.023 ±0.005 


1.022 ±0.002 


1.025 ±0.001 


1.025 ±0.0005 


2 d) 


0.140 ±0.005 


0.140 ±0.001 


0.122 ±0.015 


0.123 ±0.005 


LO: 


0.120 ±0.001 


0.1213 ±0.0004 


0.1207 ±0.0002 


0.1209 ±0.00006 



scenario 


DISASTER±± 


DISENT 


JETVIP 


MEP JET 


3 a) 
LO: 


341.2 ±1.7 
48.4 ±0.1 


339.1 ±1.2 
48.42 ±0.08 


340.1 ±0.7 
48.36 ±0.04 


331.5 ±0.4 
48.40 ±0.04 


3 b) 


as 2 a) 








3 c) 
LO: 


26.85 ±0.06 
16.94 ±0.02 


26.68 ±0.05 
16.93 ±0.02 


26.72 ±0.09 
16.93 ±0.01 


24.68 ±0.05 
16.92 ±0.01 


3 d) 
LO: 


1.998 ±0.003 
1.498 ±0.002 


1.985 ±0.003 
1.497 ±0.001 


1.995 ±0.007 
1.496 ±0.001 


1.8917 ±0.0038 
1.497 ±0.001 




scenario 


DISASTER±± 


DISENT 


JETVIP 


MEP JET 


4 a) 
LO: 


19.2 ±0.1 
11.61 ±0.04 


18.96 ±0.07 
11.57 ±0.02 


18.67 ±0.04 
11.59 ±0.01 


17.19 ±0.04 
11.587 ±0.006 


4 b) 


as 2 a) 








4 c) 
LO: 


6.42 ±0.03 
2.161 ±0.006 


6.45 ±0.02 
2.162 ±0.003 


6.40 ±0.02 
2.17 ±0.01 


6.09 ±0.03 
2.160 ±0.002 
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B Example of steering file 



The following steering file is for calculating the NLO direct cross section under HERA 
conditions in DIS with scales /i^ — Q^. 



> ========================================================================= ' 

' CONTRIBUTION' 



1 iproc [l=ep; 2=ee] 

2 isdr [1=D; 2=SR; 3=SR*; 4=DR] 
1 iborn [2->2 Born] 

1 itwo [2->2 singular contributions (for NLO)] 

1 ithree [2->3 contribution] 

isplit [gcLmma->qq term] 

iqcut [yqi-min in 2->3 matrices? l=yes; 0=no] 



INITIAL STATE' 



27.5d0 Ea [energy of lepton a] 

820. dO Eb [energy of proton/ lepton b] 

' Lepton a 

iwwa [which Weizs . -Will . : 0=ln(Q2mx/Q2mn) ; l=ln(thmax/me)] 

180. dO thmax [max angle] 

' Lepton b (only relevEint for the ee-case ) 

4.dO P2max [P2=virtuality of real photon] 

iwwb [which Weizs. -Will. : 0=ln(P2mx/P2mn) ; l=ln(thmax/me)] 

180. dO thetbmx [max angle for Weizs. -Will] 



SUBPROCESS' 



5.d0 Nf [Number of active flavours] 

0.204d0 lambda [Lambda-QCD (has to match Nf)] 

3 ialphas [QCD coupling: l=one-loop; 2=two-loop; 3=PDFLIB] 
ld-3 y-cut [phase-space-slicing parameter] 

' PDFs for the resolved contributions 

idisga [DISg -> MSbar for photon a] 

1 ipdftyp [PDF for y* (res) : l=SaS;2=GRS;3=DG;4=PDFLIB] 

2 igroupa (Param. for SaS or PDFLIB — > see meinual) 
2 iseta ( " ) 

idisgb [DISg -> MSbar for photon b] 

4 igroupb [authors of PDF on side b: y(res) or prot] 
34 isetb [Set-No] 

) Scales 

0. dO a [Scale: mu**2=a+b*C)**2+c*pt**2] 

1. dO b 
O.dO c 
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J ========================================================================= > 

' PHASE SPACE INTEGRATION' 

' ========================================================================= ' 

5d0 Q2iiiin [O.dO selects photoproduction] 

lOdO q2max 

. IdO ymin [min y] 

0.6d0 ymax [max y] 

2.25d0 ptmin [ptmin for parton in hadr cms; =Etj-min/2] 



' VEGAS and OUTPUT' 

) - — 

Ipoln [tto of Bvents producGd in phase space above] 
5 itt [no of iterations] 

'nlo.out' jfileout [Filename] 
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